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SECTION I 


INTRODUCTION 


Weather — especially precipitation — is a major consideration in the 
design of millimeter wave communications systems. Attenuation and de- 
polarization caused by rainfall, although present at lower frequencies, 
become the dominant factors affecting propagation along a ground path 
at frequencies above 12 GHz. A clearer understanding of these weather 
induced phenomena is essential so that future communications systems can 
be designed for maximum reliability and economy. Although attenuation 
has received a great deal of attention in the past 30 years, the influ- 
ence of polarization on millimeter wave propagation has been ignored 
until very recently. Practical questions which need to be considered 
relating to polarization effects are the following: (1) does a polar- 

ization exist for which the average attenuation during rainfall is a 
minimum? (2) What two orthogonal polarizations exhibit the least cross 
polarization interference? (3) Can polarization diversity be used to 
increase the reliability of a communications system by using two polar- 
izations which do not fade simultaneously? 

This thesis treats the influence of polarization on millimeter 
wave propagation from both an experimental and a theoretical viewpoint 
and provides answers to the questions listed above. It first discusses 
previous theoretical and experimental work relating to the attenuation 
and depolarization of millimeter waves by rainfall. Considerable detail 
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is included in the literature review so that the significance of the 
present work can be judged accurately. Next, a new theoretical model 
is developed to predict the cross polarization level during rainfall 
from the path average rain rate and the scattered field from a single 
raindrop. Finally, data from the VP1&SU depolarization experiment are 
presented as verification of the new model, and a comparison is made 
with other theories and experiments. 

Several features of the new theoretical model deserve mention here. 
The attenuation and cross polarization level are computed directly in 
terms of the scattered fields from the ensemble of raindrops. It does 
not depend on the difference between the attenuation of vertically and 
horizontally polarized waves (differential attenuation) to predict the 
cross polarization level during rainfall. Hence, the validity of the 
differential attenuation model is checked theoretically. In addition, 
the other theories which employ scattering and absorption cross-sections 
or effective propagation constants to predict attenuation are also 
checked. Unique aspects of the new model are: (1) spherical rather 

than plane waves are assumed, (2) the average drop diameter is used 
rather than a drop size distribution, and (3) it is simple enough so 
that the effect which changing one or more parameters has on the cross- 
polarization level is easily seen. 



SECTION II 


REVIEW OF THE LITERATURE 


2.1 Review of theoretical work 

Gustav Mie laid the foundation for theoretical work relating to the 
attenuation and depolarization of electromagnetic waves propagating 
through rain in 1908. He treated the scattering and attenuation of a 
plane wave by an ensemble of spherical particles of any material. The 
method of Mie was first applied to a rain-filled space by Ryde [2] and 
[3] during World War II, Ryde assumed that all the raindrops are spher- 
ical with the same diameter and that the rain is uniform. Medhurst [4], 
in a later review of Ryde's work, expressed attenuation as 


Attenuation = 4.343 A. 10“* dB/km 

2 it I 

where, N = number of drops per cubic centimeter 
X = wavelength in centimeters 

GO 

A^ = real part of E (a^ + P„) 


n=l 


n 


The functions a n and p n are computed in terms of the drop parameters 
and involve Bessel functions of the first kind of order n + 1/2. Although 
Medhurst [4] corrected several errors made by Ryde [3], attenuations com- 
puted using this method do not agree well with measurements. Measured 
values of attenuation tend to be well above those predicted by the theory. 
In addition to its failure to predict attenuation accurately, the 
approach used by Mie and Ryde says nothing about depolarization due to 
rainfall. Mie assumed that raindrops are spherical and spherical drops 
produce no depolarization. In order to predict depolarization it is 
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necessary to consider that real raindrops tend to be flattened in the 
vertical direction, and that they resemble oblate spheroids rather than 
spheres. 

The oblateness of the falling raindrops was considered for the first 
time by Oguchi [5] in 1960, He assumed that the eccentricity (ratio of 
the axes) of the spheroids was small and considered the oblate raindrop 
to be a perturbation of the spherical drop. He obtained the first-order 
change in the scattered field from an oblate raindrop by expanding the 
scattered field directly in a series of spherical vector wave functions 
involving powers of the eccentricity ratio v. Using the geometry shown 
in Fig. 2.1-1, Oguchi expressed the surface of the spheroid as 



where 


r = f (8) 

2 

a 



[ ; a 

^1-v Sin 2 


v < 0, prolate spheroid 
v > 0, oblate spheroid 
Then assuming that the condition 


0 


v| « 1 


held for the oblate raindrops (a condition he later found did not hold 
for the larger drops) he neglected all powers of v above the first in 
his expansion of the scattered fields. The result was a solution for 
the scattering from an oblate raindrop which gave usable results at 
least ^or small and medium sized drops. 

The next question which arises is how to extend the solution for 
a single raindrop to the actual situation of interest in which there are 
many drops in the propagation path between the transmitter and receiver. 
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r = f (e) = ■ - - a - ■ - • 

^l-v sin^ 0 
2 

where v = 1 s- 

b 

v < 0, prolate spheroid 
v > 0, oblate spheroid 

M « 1 


Figure 2.1-1. Geometry used by Oguchi in his solution for the 


scattered field from an oblate raindrop 
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Oguchi [5] and [6] used two different approaches to this problem. The 
first approach [5] was to consider the power removed from the incident 
field by a single raindrop and then sum over all the drops in a cubic 
meter in order to obtain the power removed per meter of path length. 

He defined the total cross section* Q t> of a single drop as 

_ _ total power removed from field 

"t Poynting vector of incident field 

For each raindrop Q t will consist of two parts: one part involving 

power absorbed in the drop and the other involving power scattered from 

the drop. Hence, 

^ = Qa + Q s 

and the total cross section is the sum of an absorption cross section 

and a scattering cross section. The extension to many drops is made by 

first finding the sum of the cross sections for all the drops in a cubic 

meter for a given rain rate. Denoted by <Q> this equivalent cross 

2 3 -1 

section, with units of meter per meter or meter , is defined by 

CO _ 1 

< Q > = o Q t (a) * meter 

where n(a)da is the number of drops per cubic meter which have a drop 
radius ranging between a and aH- da and Q t (a) is the corresponding cross 
section. Evaluation of the above integration requires a knowledge of 
the drop size distribution as a function of rainfall rate, which Oguchi 
obtained from Laws and Parsons [7], and a knowledge of terminal velocity 
as a function of drop radius which he obtained from Best [8], Then 
assuming that the same drop distribution is present in every cubic meter 
along the path and using the standard definition of transmission line 
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voltage attenuation coefficient, a, he derived the attenuation per meter 
as 


a = 


power lost per unit length _ <Q> meter “l 

2 • incident power at start of unit length 2 


and therefore. 


a = 4.343 x 10 3 <Q> ~ 
x km 


The second approach which Oguchi used [6] i n order to extend the 
single drop solution to the many drop problem was an application of a 
method published by van de Hulst in his book Light Scattering by Small 
Particles [9]. This method works with the scattered fields rather than 
with the power which is absorbed and scattered. An incident plane wave 
propagating in the +z direction is assumed of the form 

-j kr 

“o - S<9 -*> V* 


where S(0,$) is defined as the amplitude function for the raindrop. The 
jk factor in the denominator is arbitrary and is only inserted by van de 
Hulst for later convenience. In order to calculate the scattered field 
from an ensemble of raindrops, Oguchi assumed a plane-parallel slab con- 
taining many raindrops as shown in Fig. 2.1-2. All drops are assumed 
to be identical and identically oriented so that they all have the same 
amplitude function S(0,$). The slab has thickness 1 l and there are N 
drops per unit volume. The field at P is influenced ty scattering from 
all the drops in the slab, but the forward traveling wave is influenced 

coherently only by- those drops in the first few central Fresnel zones 
2 2 2 

(drops for which x +y << r ) . If the distance from the slab to P is 
large, then the direction is the same from each of the drops in the 
central Fresnel zones to P and the amplitude function, 8(0,$) = S(0,0), 
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. Oguchi’s application of van de Hulst's method to find 
the attenuation caused by many raindrops. 
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is the same for all of these drops. The total amplitude at P is then 

u = u o {1 + S(0,0) j e~ ik (* + y )/2r dx dy dz) 

where the integration is over the entire slab, but where it is assumed 

2 2 2 

the S(0,0) is zero for drops for which the condition x +y << r does 
not hold (see van de Hulst [9] page 34). This is a form of Fresnel's 
integral and can be integrated to give 

u - u {i - ~ m s(o,o)} 

o K z 

This result may be formally represented as the influence of a complex 
refractive index of the medium as a whole. If the slab is replaced by 
a slab of homogeneous material with complex refractive index m, m close 
to 1, the amplitude of the wave is changed by the slab in the proportion 

e -jk£(m-l) = 1 _ j k£ ( m _i) 

Comparing the above equation with the equation for u it is seen that 

m = 1 - j S(0) 2ir N k~ 3 

Separating m into a real part, n, and an imaginary part, n', we have 

n = 1 + 2 tt k" 3 I <S(0)} 
m 

n' = 2tt k 3 Re {S(0)} 

The real part, n, determines a phase shift of the wave traveling through 
the rain while the imaginary part, n', determines the attenuation. 
Oguchi, and more recently Setzer [10], have used this method to compute 
the attenuation and phase shift of a wave traveling through a rain- 
filled space. Oguchi [6] has computed the amplitude function S(0,<P) 
for both vertical and horizontal polarization and shown that a 
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horizontally polarized wave will be attenuated more by rain than will a 
vertically polarized wave. This is a direct consequence of the oblate 
shape of the raindrops. 

This difference in attenuation between vertically and horizontally 
polarized waves which Oguchi calculated has been the basis for all 
theoretical work on rain depolarization to date. Application of Oguchi' s 
results to the depolarization problem was first made by Thomas [11] at 
Bell Labs, who showed that a difference in attenuation for waves polarized 
along the major and minor axes of an oblate raindrop would lead to de- 
polarization of any wave not polarized parallel to one of these prin- 
cipal axes. The major and minor axes of the drop are the greater and 
lesser axes of symmetry of the vertical cross section of the drop, 
respectively. The depolarization effect is illustrated in Fig. 2.1-3 
where E T represents a transmitted wave whose polarization is not parallel 
to either principal axis. E^ can be expressed as a component E^ along 
the major axis and a component E 2 along the minor axis. As E-^ and E 2 
propagate, they will be attenuated but their direction will not be 
changed since they are parallel to the principal axes. and E^ re- 
present E^ and E^ at the time they reach the receiver. Since E^ has 
been attenuated more than E^, E^ will be proportionally smaller than E^ 
and the resultant will not be parallel to Ej. Hence E^, has been 
depolarized or rotated in the direction of the minor axis of the rain- 
drop. With reference to Fig. 2.1-3 the cross polarization level (XPOL) 
is computed as follows. 


Define, Ay «* attenuation of vertically polarized wave 

A^ <= attenuation of horizontally polarized wave 



1 


r 



differential attenuation 
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n = angle between Incident E field and vertical 

E 2 = \ l S T l cos 1 

E { = Ar l\l sin ^ 


QXU 1 | A J - 


e r = (A v I\1 cosn y + Ajj |e t | 

= Ay | E t | COS 2 n + An 1^1 sin 2 n 

= (Ay I Ej, I cos n y + Ajj I ^ | sin n x) • (sin n y - COS n x) 
= Ay I E t I cos n sin n - A^ |e t | cos n sin n 

E^ Ay cos n sin . n - A^ cos n sin n 

__ ■ - 2 “T ~~2 

EL Ay cos n + Ajj sin n 
cos n sin n . (1 — — ) 


2 . lA H . 2 

cos n + -r— sin n 

Ajj 

cos n sin rj (1 — -r~) 
\ 

1 + sin n (~ - 1) 

*V 


e r 

XPOL = 20 log 1(J ~ - 20 log 1Q 
• \ 


Ajj 

cos n sin n (1 j— ) 

\ 


2 \ 

1 + sin n (t 1) 


These equations represent the state-of-the-art for computation of the 
cross polarization level during rainfall. As is now evident, most of 
the theoretical work has considered attenuation produced by rain rather 
than depolarization. The differential attenuation model is the only 
work done to date on the theoretical problem of depolarization and many 
authors have used it to make quantitative estimates of cross polarization 
levels. The position of the differential attenuation model is that of 
any new theory which solves a previously unsolved problem — it must be 
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proved or disproved by further work. In addition to checking the 
validity of the differential attenuation model, the new model developed 
in this thesis overcomes several limitations of the old theory. The 
most notable contribution is that the cross polarization level is com- 
puted directly in terms of the summation of the scattered fields from 
the ensemble of raindrops rather than being computed indirectly from 
attenuation predictions. This makes it much easier to see how changing 
various parameters such as drop size, path length, drop orientation, 
rain rate, terminal velocity, etc., affect the cross polarization level. 
The use of a drop size distribution is found to be unnecessary. The 
uncertainty as to the number of oblate raindrops severely limits 
accurate prediction of cross polarization and precludes the use of 
anything other than an average drop size. 

The ultimate test of a theoretical model is whether or not it 
accurately explains experimental observations. In order to correlate 
theory and experiment, it is imperative that the experimental data be 
taken under well defined conditions of rainfall rate, drop orientation, 
etc. The next section discusses the previous experimental work 
relating to depolarization and it is concluded that no experiment has 
been performed which serves as more than a crude check of a theoretical 
model. These experiments must be reviewed in order to point out their 
limitations and to decide how a depolarization experiment should 
performed. » * 

2.2 ' Review of experimental work 

In 1965, Medhurst [4] published an excellent review of the available 
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results of experimental attenuation measurements using millimeter waves. 

He compared the theoretical attenuations computed by Ryde [2] and [3] 
to experimental data concerning attenuation by rainfall. His conclusion 
was that "there is a marked tendency for observed attenuations to fall 
well above levels which, according to the theory, cannot be exceeded." 

Some of the pre-1965 experiments merit review here in order to learn how 
such experiments should be conducted. 

Typical of the first experiments to measure attenuation due to rain- 
fall was one performed in 1942 by Adam [12], He used a frequency of 
31.25 GHz and a path length of 2 km. He had only one rain gauge along 
the entire path and it was located at the transmitter end. Semplak [24] 
found that rain rate can vary as much as 100 mm/hr in a 0.5 km distance 
and that rain gauge spacing should be less than 0.5 km for accurate rain 
rate determination. This means that the path average rain rate cannot 
be inferred from Adam's data making his attenuation measurements of little 
value. Other early experimenters who used too few rain gauges were 
Hathaway [13], Rado [14], Okamura [16], and Funakawa [17], and their 
measured attenuations are equally questionable. 

The first experiment utilizing a sufficient number of rain gauges 
to accurately measure the rainfall rate along the path was conducted in 
1947 in Hawaii by Anderson [18], The frequency used was 24 GHz and the 
path length was 1.95 km. Nine rain gauges were distributed along the 
path. Medhurst [4] considers this to be the most carefully performed 
of any experiment which he reviewed. Figure 2.2-1 shows the data taken 
by Anderson and it can be seen that the measured attenuation is well 
above the theoretical curves as calculated by Ryde [3], This may be 
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Figure 2,2-1. Measured rainfall attenuation from Anderson [4] 
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due in part to the tropical climate which produces orographic rain 
rather than convective thunderstorms. 

Perhaps the most advanced measurements of attenuation prior to 1965 
were made in 1959 by Okamura [15] at the Radio Research Labs in Tokyo. 

He used a frequency of 34.884 Gllz and a path length of 400 meters. How- 
ever, a reflection method was used so that the radio wave traversed the 
path twice making an equivalent path length of 800 meters. One rain 
gauge was used which was read every 15, 30, or 60 seconds depending on 
the rainfall rate. The unique feature of Okamura' s experiment was that 
he used both vertical and horizontal polarization, transmitting both 
during a five minute interval. Figure 2.2-2 shows the attenuation data 
which he obtained for each polarization. The theoretical limits are 
from Ryde [6], As can be seen the attenuation is greater for horizontal 
polarization then it is for vertical polarization and that both tend to 
be at or above the maximum possible attenuation as predicted by theory. 

In 1967 Semplak [19] began an extensive set of rainfall attenuation 
measurements covering a period of four years. The first experiment was 
designed to measure the attenuation of a vertically polarized wave using 
a frequency of 18.5 GHz and a path length of 6.4 km. Figure 2.2-3 shows 
the attenuation data which Semplak obtained for a total of 21 storms. 

The rain gauge spacing of 1.6 km is a possible explanation for the large 
amount of scatter in the data. 

In 1969 Semplak [20] measured the difference^ between the. rainfall 
attenuation of vertically and horizontally polarized waves using a 
frequency of 30.9 GHz and a path length of 1,89 km. He transmitted 
a field polarized at an angle of 45° from the vertical and switched the 
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Rainfall rate in mm/hr 
(a) Vertical Polarization 



Rainfall rate in mm/hr 
(b) Horizontal“Polarizatibn 

Figure 2,2-2. Measured rainfall attenuation from Okamura [4] 
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Rainfall rate in mm/hr 


Figure 2,2-3. Measured rainfall attenuation from Semplak 

for a 6.4 km path at 18.5 CHz [19] 

(Copyright 1969, American Telephone and Telegraph 
Company, reprinted by permission) 
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receiving antenna between horizontal and vertical polarization at a 
rate of 329 Hz using a ferrite switch. He recorded the sura and dif- 
ference of the horizontal and vertical fields and thus obtained the 
differential attenuation as a function of total path attenuation. His 
cumulative data plot is shown in Figure 2.2-4 with horizontal polar- 
ization being attenuated more than vertical. 

In 1969 Skerjanec and Samson [22] made simultaneous measurements 
of rainfall attenuation at 10 and 14.43 GHz. They had 20 tipping-bucket 
rain gauges on a 4.7 km path. Measured attenuation (in dB) at 14.43 GHz 
was approximately twice that observed at 10 GHz when the rainfall rate 
exceeded 50 mm/hr. The maximum attenuation observed was 24 dB at 10 GHz 
and 48 dB at 14.43 GHz corresponding to a rain rate of 137 mm/hr. 

Readings were recorded only when the rain rate was relatively uniform 
at all gauges over a selected interval (1-5 min) . 

In 1970 Semplak [21] made simultaneous measurements of rainfall 
attenuation at 18.5 and 30.9 GHz. By examining the ratio of the 30.9 GHz 
attenuation to the 18,5 GHz attenuation during rain, he concluded that 
the Law and Parsons [7] drop size distribution is adequate to predict 
the rainfall attenuation observed with terrestrial microwave communications 
systems. 

Since May 1971, Watson [23] at the University of Bradford in England, 
has been making measurements of attenuation and cross polarization due 
to rain and snow at a frequency of 11 GHz. He is using ten rain gauges 
along a 13.5 km path. Initially he transmitted only vertical polarization 
and employed co-polar and cross-polar antennas at the receiver. Recently 
a separate horizontally polarized antenna was added at the transmitter 
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Figure 2.2-4. Measured differential attenuation from Semplak 
for a frequency of 30.9 GHz and a path length 
of 1.89 km. [20] 

(Copyright 1970, American Geophysical Union, 
reprinted by permission) 
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and the polarization transmitted is being switched between vertical and 
horizontal at 10 Hz. Figure 2.2-5, which is a plot of cross polariza- 
tion vs. attenuation, is typical of the data which Watson has taken. 

As can be seen, there is no correlation between the two and, since rain- 
fall rate and attenuation are highly correlated, there can be no cor- 
relation between cross polarization and rainfall rate. This is due, 
as Watson says, to the uncertainty as to the canting angle of the rain- 
drops. Theoretical curves of cross polarization vs. attenuation for 
canting angles of 5 and 45 degrees are shown in Figure 2.2-5. The problem 
is that, for a canting angle of zero degrees, the cross polarization 
isolation is infinite regardless of the attenuation if the transmitted 
wave is either vertically or horizontally polarized. Hence vertical or 
horizontal polarization is extremely sensitive to the raindrop canting 
angle (which is never known). This makes it very difficult to accurately 
check a theoretical model of cross polarization due to rainfall with data 
obtained from an experiment using vertical and/or horizontal polarization. 
The new theoretical model, which will be developed next, predicts the 
optimum polarization to use for the best correlation of theory to 


experimental data. 
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Figure 2.2-5, Measured rainfall attenuation vs. cross 
polarization from Watson. 



SECTION III 


DEVELOPMENT OF THE THEORETICAL MODEL 

3.1 Assumptions 

Assumptions made in the theoretical model development are: (1) a 
uniform rain along the path, (2) that all raindrops, regardless of their 
shape, have the same volume equal to the volume of the most frequently 
occurring drop, (3) a rain-filled space which is large-scale homogeneous, 
(4) transmitting and receiving antennas with sufficient beamwidth so 
that small variations in beamwidth do not affect the received signal 
(this implies that many Fresnel zones are included in the common volume 
of the beams), (5) both the incident and scattered fields are spherical 
waves, and (6) the scattering function for a single drop is known. 

3.2 An Essential Theorem 

Translocation theorem: For a rain cell of width Afi, occurring along 

a propagation path, the total field at the receiver is independent of 
the actual position of the cell along the path and depends only on the 
cell width. 

The above theorem is central to the theoretical model development 
and its validity will be demonstrated first. The technique to be 
employed is the method of Fresnel zone analysis which has been used 
extensively in the solution of diffraction problems, Livingston [25] 
gives a detailed justification for this method which will be reviewed 
briefly here. 
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Consider an arbitrary plane perpendicular to the propagation path 

in a homogeneous medium as shown in Figure 3.2-1. Point 0 is the inter- 

cept of the plane with a straight line connecting the antennas. By 
Huygen's principle the field at the receiver may be obtained by con- 
sidering each point in the plane to be the source of a new spherical 

wave with an amplitude and phase equal to the amplitude and phase of the 

transmitted field at that point. The total field at the receiver is 
then the sum of the spherical waves from each point in the plane. The 
wave from point 0 will arrive at the receiver first since it travels the 
shortest distance. The phase of the waves from all other points will 
lag behind that of the wave from point 0, and this phase lag will increase 
as points farther and farther from 0 are considered. The locus of points 
in the plane for which the phase lag is exactly tt radians determines a 
circle with center at 0 and radius F^. This circle defines the boundary 
of the first Fresnel zone. All points at a distance less than F^ from 
point 0 are in the first Fresnel zone. Similarly, circles with radii 
* 2 * ^ 3 » etc * are defined for phase shifts of exactly 2 tt, 3tt, etc. with 
the area between F^ and defined as the second Fresnel zone, that 
between ¥ ^ and F^ the third Fresnel zone, etc. The total field at the 
receiver may be expressed as a summation of the fields from each of the 
Fresnel zones. Since there is an infinite number of Fresnel zones, an 
infinite series must be summed. Since the phase of the wave is changing 
by -ir radians as the summation proceeds from one zone to the next, the 
signs of the terms of the infinite series alternate. Livingston shows 
that each of the Fresnel zones has the same area and therefore the same 



Arbitrarily located plane 
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Figure 3.2-1. Definition of Fresnel zones 
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number of points to contribute to the field, but the contribution of the 
outer zones is increasingly smaller due to the greater distance and cor- 
respondingly greater free-space loss. The summation of the fields from 
each of the Fresnel zones is, therefore, the summation of an infinite 
geometric alternating series. Since the sum of such a series approaches 
half the value of the first term of the series when the value of its 
common ratio approaches unity, the total field at the receiver may be 
found by calculating the contribution of only those points in the first 
Fresnel zone and taking one-half of the result. 

Now for a proof of the theorem, consider the geometry as shown in 
Figure 3.2-2 where the first Fresnel zone is defined as the locus of 
points for which the phase lag is it radians. The total path length 
between the transmitter and receiver is L, and A £ is an elemental length 
of the path at a distance £ from the transmitter. AV is the elemental 
volume defined by A£ and the boundary of the first Fresnel zone. There 
is assumed to be rain in the volume AV with no rain along the rest of 
the path. The radius of the first Fresnel zone at £ is given by 
Livingston [25] 

r . <» - A!/' 2 . 

The vertical axis in Figure 3.2-2 is greatly expanded relative to the 
horizontal axis. At 17.65 GHz the radius r is less than one meter at 
midpath for a path length of _1. 43 km. ^Assuming. spherical waves- are 
transmitted, the incident field at the receiver will be inversely pro- 
portional to the path length L while the incident field at £ will be 
inversely proportional to £. 
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Figure 3.2-2. Geometry Used in Proof of Theorem. 
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The scattered field at the receiver from the raindrops in the volume AV 
will be directly proportional to the incident field at S and directly 
proportional to the total number of raindrops in AV which in turn is 
proportional to AV. The scattered field at the receiver will be in- 
versely proportional to the distance L -S from AV to the receiver. 
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Hence the scattered to incident ratio is independent of the position of 
a rain cell along the path and is independent of the length of the path. 
It depends only on the width of the rain cell. If rain is falling along 
the entire path, then AS=L and the ratio depends on the path length. 


3.3 . Model iDevelopment • - - - - - - 

3.3.1 Path Division 

Consider the rain-filled volume bounded by the first Fresnel 
zone as shown in Figure 3.3.1-l(a). Divide the total number of drops 
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Figure 3. 3.1-1. Division of Path When 
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into N groups by dividing the path into N equal segments of length 
L/N. Place the drops in each group in a plane at the center of each 
segment as illustrated for N=4 in Figure 3.3.1-l(b). It follows 
from the translocation theorem of Section 3.2 that the contribution 
to the scattered field at the receiver is the same for each plane 
of drops and is independent of the actual position of the plane along 
the path. However, it must be noted that planes near the center of 
the path contain more drops than those near the ends since the volume 
of a segment near the center is larger than the volume of a segment 
near either end. Moving a plane of drops from one position to 
another requires changing the number of drops by the ratio of the 
volumes of segments at the two positions in order to maintain the 
same contribution to the scattered field at the receiver. 

Next assume that all planes are moved to new positions near 
the center of the path as shown in Figure 3.3.1-l(c). All planes 
will be at the same distance L/2 from the receiver and all planes 
must make equal contributions to the scattered field (since they 
all made equal contributions in their former locations). It follows 
that all planes must have the same number of drops and are identical, 
Therefore one need only solve for the scattered field from a single 
plane. 

3.3.2 Scattering from a Single Plane 

The following complex forward scattering coefficients are 
defined for a single plane of drops in terms of an incident (main) 
polarization (1) and an orthogonal (cross) polarization (2) . These 
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coefficients will be the same for all planes. 


'll 


scattered E field with polarization (1) produced by an 
incident E field with polarization (1). 
scattered E field with polarization (2) produced by an 
incident E field with polarization (1) . 
scattered E field with polarization (1) produced by an 
incident E field with polarization (2). 

S 22 scattered E field with polarization (2) produced by an 
incident E field with polarization (2). 

The scattered field at the receiver from a single plane of drops 
can be found by summing the scattered fields from each of the drops 
in the plane. With reference to Section 3.2 it is necessary to con- 
sider only those drops in the first Fresnel zone and to multiply 
the result by one-half. Oguchi [26] has computed complex forward 
scattering coefficients for a single raindrop. In terms of for- 
ward scattering functions (fpq) computed from Oguchi' s functions, 
(see the Appendix) the scattering coefficients for a plane of drops 
will be of the form 

-j k r 


S = r H 1 f - 

pq 2 p pq r 

where E^ is the incident field with polarization p, 

r is the distance from the drop to the receiver, 

f is the ratio of the scattered E field 4*ith polar- 

pq 

ization (q) to the incident E field with polar- 
ization (p) for a single drop, 
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and where the summation is to be performed over those drops in the 
first Fresnel zone. The scattering function f is assumed constant 
for all drops in the plane which have the same shape. As discussed 
in Section 4, the total drop population is assumed to consist of 
40% identical oblate drops and 60% identical spherical drops. A 
constant f for all drops of the same shape requires that all 
oblate drops be identically oriented and that the radius of the 
first Fresnel zone be very much less than the path length so that 
only forward scatter need be considered. The r in the denominator 
of the previous equation may be replaced with L/2 yielding 

E 1 f 

S --JL-BS. e e “ jkr 

pq L 


Assuming a uniform distribution of dropB within the first Fresnel 

zone, the above equation is easily evaluated. The result of the 

2 

summation is an imaginary constant, K = -j — , times the number of 
drops, D, in the plane (see the Appendix). Hence, 



E 1 f 

jb— p.a 
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J ttL p 


f 

pq 


3.3.3 Effect of H Planes 

The effect of N planes of drops in series may be found by 

considering the incident field to interact with the first plane 

producing a new field which then interacts with the second plane, 

etc. until all N planes have been considered. This may be expressed 

in closed form in terms of the scattering coefficients (S 's). It 

pq 

is assumed that the incident field at midpath has amplitude two so 
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that the incident field at the receiver has unity amplitude. It 
is also assumed that the contribution from the S 2 ^ coefficient 
(scatter from cross to main polarization) is small and that S 2 2 
is zero. This is valid for cross polarization levels less than 
about -5 dB. 

As shown in Figure 3. 3. 3-1, the total fields at the receiver 
are obtained by induction from the effect of a single plane. After 
interaction with the first plane of drops the total main polar- 
ization field will consist of the sum of the incident field at 
the receiver (unity) and the scattered field S^. The total cross 
polarization field E^ will consist only of the scattered field S 12 
since the incident field in this polarization is zero. The fields 
leaving plane 1 then interact with plane 2. The main polarization 
component will be multiplied by another (1 + S.^) factor. The cross 
polarization component will consist of two parts: (1) the previous 

cross polarized component E^ which will be attenuated by a (1 + S 22 ) 
factor, and (2) a new contribution from the previous E^ equal to 
times S^ 2 « This method is continued for all N planes with the 
following result as shown in Figure 3. 3. 3-1. 

" (1 + S l/ 

E X = S 12 t(1+S 22 )N ~ 1 + (1+S 22 ) N ~ 2 ( 1 +S ii> + •••• 

' + (1+S 22 ) (1+S 11 ) N “ 2 + (l+Sj^)^ 1 ] 

The above described method forms the basis for the computer 
program which is listed in the Appendix. This program was used to 
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Figure 3. 3.3-1. Effect-of-N planes of drops in series. 



make calculations of expected cross polarization and attenuation 
levels during rainfall which are compared to experimental data 
from the VPI&SU depolarization experiment in the next section. 

In Section 5 a comparison is made with other theories and experiments. 



SECTION IV 


EXPERIMENTAL VERIFICATION 

4.1 VPI&SU Depolarization Experiment 
4.1.1 General System Description 

The depolarization experiment sponsored by NASA and con- 
ducted at Virginia Polytechnic Institute and State University has 
been described in detail elsewhere [27] and will be reviewed only 
briefly here. The experiment is designed to measure depolarization 
and attenuation effects of rainfall on millimeter wave propagation. 
The path is 1.43 km in length and is located on the VPI&SU campus 
in Blacksburg, Virginia. The measurements are being made at a 
frequency of 17.65 GHz with two separate transmit and receive 
channels orthogonally polarized at i45° from the vertical. 

Figure 4. 1.1-1 presents a block diagram of the overall 
experiment. It consists of (1) an RF system (millimeter wave 
transmitter and receiver), (2) transmitting and receiving antennas 
(dual polarized), (3) a weather monitoring system with rain gauges 
and wind sensors, (4) a digital control, processing, and data 
storage system built around a Raytheon PB 440 computer, and (5) a 
link via magnetic tape to an IBM 370/155 computing system for off- 
line data processing. 

An important feature of this experiment is the absence of 
possible multipath effects which could produce cross polarization. 
The common volume of the main beams of the antennas does not 
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intercept the ground or any other obstacle. The angle to the 
ground at midpath from either antenna is 2° and the angle from the 
mainbeam maximum to the first null of the radiation pattern is 
about 1°. Therefore only sidelobes intercept the ground and any 
multipath effects are more than 40 dB below the direct signal. 

4.1.2 Data Processing 

4. 1.2.1 On-Line Processing 

A Raytheon PB 440 computer assisted by a special- 
purpose controller operates the experiment, acquires data, 
and does preliminary data processing. The experimental 
control program maintains the system in the proper operating 
mode for current weather conditions and signal behavior. 

The clear weather operating mode is called mode 0, and in 
it the +45° transmitter channel operates continuously while 
the computer monitors the +45° to -45° cross polarization 
level and the +45° direct attenuation. Both receiver 
channels are sampled at 10 second intervals while wind 
velocity and transmitter power are sampled every 100 seconds. 
If the cross polarization level (in dB) changes by more 
than 2% or if one of the rain gauges reports precipitation, 
the system begins operating in mode 1. During mode 1 
operation, transmission is sequenced at 4 second intervals 
from the +45° to the -45° channel and then to both channels. 
Receiver sampling occurs at 1 second intervals and wind 
velocity is sampled every 4 seconds. Mode 1 operation 
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continues until the precipitation rate falls below 6 mm/hr 
or until the cross polarization level stabilizes. At this 
time, mode 2 operation is begun with transmitter switching 
at 10 second intervals and receiver and wind sampling at 
2 and 10 second intervals respectively. Mode 2 operation 
continues until the precipitation rate falls below 3 mm/hr. 

The system then enters mode 3 with transmitter switching at 
100 second intervals and receiver and wind sampling at 10 
and 100 second intervals respectively. When the precipitation 
rate falls below 2 mm/hr, the system re-enters mode 0 
operation. In all modes there is a low pass filter (time 
constant = 0.4 seconds) at the input to the A-D converter 
which surpresses fast scintillations of the signals and 
insures that average values are sampled. 

When a new data point enters the computer, a pro- 
gram locates the last two values stored for that input. 

If the new value and the last value differ by more than 1% 
the new value is stored. If the difference between the new 
value and the last value is less than 1%, the new value is 
compared to the next to last value. If these differ by more 
than 1%, the new value is stored; if this difference is less 
than 1%, the last value is discarded and the new value takes 
its place. 

4.1. 2. 2 Off-Line Processing 


An IBM 370/155 computer program has been developed 
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which processes, analyzes, and plots the accumulated data 
from any number of selected storms. These data are rain 
rates from each gauge plus quasi- instantaneous (i.£. short 
integration time) samples of the analog signal levels during 
a storm. The latter are stored at essentially regular times 
while the intervals between successive rain gauge trips are 
random. Before data from different inputs can be compared 
the computer must generate a time-function representation 
for each variable. These time functions are then averaged 
over selected time intervals to generate the average signal 
levels, rain rates, etc., required by steady state theory. 

4.2 Comparison of Theory and Experiment 

The data from the VPI&SU depolarization experiment have been des- 
cribed in detail elsewhere [28] . The intent of this section is to pro- 
vide a review of the data sufficient to check the validity of the de- 
polarization model developed in Section 3. Experimental results are 
presented using graphs which also contain theoretical predictions based 
on the model of Section 3 and calculated using the computer program 
listed in the Appendix. These figures make the comparison between theory 
and experiment evident at a glance. The main test of the theory is 
the accurate prediction of the cross polarization level as a function 
of rainfall rate and most of the figures are of this type. However, 
plots of attenuation versus rainfall rate and cross polarization versus 
attenuation are included since they provide additional support for the 
model. For all of the data presented here a 15 second time average 
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(performed off-line) was used for both the rain rate and the received 
signal levels. In addition, all data points for each integer value of 
rain rate were averaged to produce a single data point for that rain 
rate. Standard deviations were also computed for each data point. The 
range of standard deviations is stated in the text for each figure but 
it is not plotted since it tends to obscure the data points. 

The first storm after the experiment became operational occurred 
on August 4, 1972. Figure 4.2-1 shows the cross polarization level for 
this storm. There is considerable scatter in the data but the theoretical 
curve is certainly representative of the average level observed. The 
scatter is most likely due to an inaccurate measurement of the rainfall 
rate since only two raingauges were in operation at the time. As dis- 
cussed in Section 2.2, an intergauge spacing of 1.43 km is not sufficient 
for an accurate measurement of rainfall rate. The standard deviation of 
the data points in Figure 4.2-1 is about 5 dB and all but two points are 
closer than this to the theoretical curve. 

Figure 4.2-2 presents data taken on August 17, 1972. This storm 
provided most of the data which were taken during 1972. The peak rain 
rate of 105 mm/hr was the maximum rain rate observed. As can be seen, 
there is good agreement with the theoretically predicted cross polar- 
ization level. Since there were still only two raingauges along the 
path, scatter is to be expected in this data also. The standard devi- 
ation for these points is" a little less than 5 dB. 

Data from the storm of September 29, 1972 is shown in Figure 4.2-3. 
There is excellent agreement with theory except for rainrates less than 
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20 mm/hr. This is because the residual cross polarization level of the 
system at that time was about -30 dB. These points do not represent an 
effect of rainfall along the path. The scatter is greatly reduced in 
this data since all five rain gauges were in operation by this time. 

The standard deviation is 6-8 dB. 

Figure 4.2-4 shows data taken on October 27, 1972. There is good 
agreement with theory and a lack of scatter which again may be attributed 
to the operation of all five rain gauges. Standard deviations were not 
computed for these points. 

Figure 4.2-5 presents the average cross polarization level for the 
six heaviest storms occurring between August 1 and December 31, 1972. 

Also shown are theoretical curves for both 40% and 100% oblate rain- 
drops. The number of oblate raindrops assumed is seen to have a decided 
effect on the predicted cross polarization level. Previous to the 
present work, it has been assumed when computing depolarization that 
100% of the drops have the same oblate shape, although Jones [29] has 
found that a real rain consists of drops with a variety of shapes. He 
found the following distribution of shapes for all sizes (greater than 
1.9 mm equivalent spherical diameter) taken together: 


spherical 

32% 

oblate 

28% 

prolate 

18.5% 

irregular 

21.5% 


All drop shapes except spherical could contribute to the cross polar- 
ization level. About one-half of the prolate drops would be expected 
to be aligned with their elliptical cross section facing the incident 




27, 1972 
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field and therefore produce depolarization. This would increase the 
effective percentage of oblate drops to about 37%. Since a small but 
unknown contribution would be made by the irregular drops, the effective 
percentage of oblate drops was fixed at 40% for this analysis. For 
percentages near 40% an increase of 5% in the assumed number of oblate 
drops increases the cross polarization level by about 1 dB for all rain- 
rates. From Figure 4.2-5, it is evident that an assumption of 40% 
oblate drops provides good correlation between the theoretical model 
and the experimental observations. 

Another way to compare the theory with the experimental data is 
to plot attenuation versus rain rate and this is done in Figure 4.2-6 
for the storm of August 17. This figure contains the same type of 
scatter as the figures illustrating the cross polarization levels 
although the theoretical curve represents the body of the data very 
well. This scatter is due to an inaccurate determination of the rain 
rate (due to nonuniform rain along the path) and can be shown by 
plotting attenuation versus cross polarization thereby eliminating 
rain rate as a parameter. The result (for the storm of August 17) 
is shown in Figure 4.2-7 where the scatter is no longer present and 
an excellent agreement is obtained between theoretical prediction and 
experimental observation. 

From the various figures which have been presented here, it is 
clear that the new theoretical model of rain depolarization developed 
in this thesis accurately predicts the average cross polarization and 
attenuation levels which were observed in the VP1&SU depolarization 
experiment. In Section 5 a comparison is made with other theoretical 
and experimental results. 





SECTION V 


RESULTS AND CONCLUSIONS 

5,1 Comparison with Previous Theories 

This thesis has demonstrated that the method of Fresnel zone 
analysis can be applied with good results to the problem of predicting 
the cross polarization level during rainfall. Since this approach 
differs fundamentally from that taken by previous investigators, it 
provides, for the first time, a theoretical means for checking their 
results. The most noted work to date is by Oguchi [26] who has cal- 
culated attenuation and phase rotation for both vertically and horizon- 
tally polarized waves as a function of rain rate. A comparison with 
his predicted attenuation at 19.3 GHz for a 1.43 km path as a function 
of rain rate is shown in Figure 5.1-1, As can be seen, the maximum 
difference between the two curves is about 0.5 dB at 150 mm/hr. Figure 
5.1-2 compares Oguchi 's computed phase rotations at 19.3 GHz for a 
1.43 km path and vertical polarization with those obtained using the 
new model. The agreement here is not so good as in Figure 5.1-1 but 
there is as yet no experimental data to support either curve. Another 
theory which can be checked in light of the present work is the dif- 
ferential attenuation model of depolarization. This method was developed 
by Thomas [11] and until now was the only way in which the cross polar- 
ization level during- rainfall-could be-predicted. A comparison with 
the differential attenuation prediction using Oguchi' s values of 
attenuation is shown in Figure 5.1-3. For the new model, curves are 
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Figure 5.1-1. Theoretical Prediction of Rainfall Attenuation at 19.3 GHz 
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5.1-2. Theoretical Phase Rotation for Vertical Polarization at 19.3 GHz 
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shown for both 40% and 100% oblate raindrops. As discussed in Sec- 
tion 4.2, the actual effective percentage of oblate drops is near 40% 
and this is the value which should be used in computations . The 40% 
curve for the new model is in good agreement with experimental data. 
However, when using the differential attenuation method, 100% oblate 
drops is assumed and the middle solid curve of Figure 5.1-3 is obtained. 
This curve agrees reasonably well with experiment, predicting values 
that are only slightly high. But the difference in the predicted cross 
polarization level when as sinning 40% and 100% oblate drops is about 
8 dB. If the differential attenuation curve is shifted down 8 dB, the 
dashed curve of Figure 5.1-3 results which is well below experimental 
observation. Hence, it is concluded that the differential attenuation 
model of depolarization does not agree well with experimental measure- 
ments if a realistic percentage of oblate drops is assumed. 

Another conclusion which can be drawn is that the use of a dis- 
tribution of raindrop sizes is unnecessary in order to get good agree- 
ment between theory and experiment. A far more important factor that 
must be considered (which has been neglected previously) is the 
effective percentage of oblate raindrops. 

5.2 Depolarization of Arbitrarily Polarized Waves 

The new model of depolarization developed in this thesis may be 
coupled with the general theory of depolarization of Beckmann [30] to 
obtain the depolarization by rain of an arbitrary elliptically polarized 
wave, Beckmann’s theory states that the depolarization of an arbitrarily 
polarized wave can be found if the depolarization of two orthogonal 
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linearly polarized waves is known. His method is illustrated in Fig- 
ure 5.2-1 where E^, represents the electric field vector of a transmitted 
arbitrarily polarized wave at a particular instant of time. E^ can 
be resolved into components E Tx and E^,^ along the x and y axes respec- 
tively. As E^ x propagates toward the receiver, a part of its energy 
will be converted to a y-polarized wave which is called The 

remainder of the energy of E^ which reaches the receiver will still 
be x-polarized and is denoted by Ej^. Similarly, E^ produces a cross 
polarized component E^ x and a direct component E^ . Hence each com- 
ponent of the received wave is composed of the sum of a direct wave and 
a cross polarized wave. In equation form. 


EL = E„ + E_ 

Rx Rxx Ryx 


and \y ~ \yy + \ 


yx 


The following depolarization coefficients are defined: 
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The r’s are complex quantities and require that both the magnitudes and 
phases of the received signals be known. The new model provides the 
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magnitude and phase for both the direct and the cross polarized wave and 
allows a computation of Beckmann’s depolarization coefficients for a 
wave propagating through rain. Beckmann shows that 


r p. + r 

P = _yy_ 1 y* 

2 r + r p ■ 

xx xy 1 

where P^ and P^ are the complex polarization factors of the transmitted 
and received waves respectively. The received polarization during rain 
for any transmitted polarization can be computed using the above equation 
and T's calculated for the new model. 


5.3 Preferred Polarization for Depolarization Measurements 

Vertical and/or horizontal polarizations have been the choices of 
most experimenters measuring depolarization by rain. At first thought, 
this seems to be the logical choice, but it leads to serious problems 
when it is desired to correlate the experimental data with a theoretical 
model. The reason is the uncertainty as to the canting angle of the 
raindrops. The canting angle is the angle which the minor axis of an 
oblate raindrop makes with the vertical and is zero if the rain is 
falling straight down. But rain usually does not fall straight down 
and typical canting angles are in the range from -15° to +15°. Figure 
5.3-1 shows the cross polarization level for two paths as a function 
of the tilt angle, where the tilt angle is defined as the angle between 
the incident E field vector and the minor axiB of an oblate raindrop. 

The cross polarization “level” for' tilt angles of 0° or 90° is -«*. If the 
transmitted field is polarized either vertically or horizontally, then 
the canting angle is equal to the tilt angle and a rain with zero canting 
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angle will produce no depolarization (since the cross polarization level 
is However, small variations in the canting angle about 0° or 90° 

produce large changes in the predicted cross polarization level. Since 
the magnitude of the canting angle could be anywhere from 0° to 15° or 
more, all that can be done when vertical or horizontal polarization is 
used is to place broad limits on the expected cross polarization (with 
being the lower limit). 

In contrast, the curves of Figure 5.3-1 are relatively flat for 
tilt angles of 45° to 60°. Consider the lower curve which is for the 
1.43 km path used in the VPI&SU depolarization experiment. If a trans- 
mitted field polarized at 45° from the vertical is assumed, then the 
tilt angle is 45° when the canting angle is zero. For canting angles 
of ±15° the tilt angle will range from 30° to 60°. The predicted cross 
polarization level ranges from a maximum of -13 dB at 55° to a minimum 
of -16 dB at 30°. The variation is only 3 dB for the expected range 
of canting angles. The advantage of using +45° polarization is evident 
after a comparison of Figure 2.2-5 on page 22 to Figure 4.2-7 on page 50. 
The former shows the attenuation and cross polarization data that Watson 
has taken using vertical and horizontal polarization. As can be seen, 
there is no correlation between theory and experiment. The latter shows 
data from the VPI&SU depolarization experiment using i45° polarization. 
Excellent correlation is obtained in this case. 

It is apparent that vertical and horizontal polarization are the 
poorest choices for a depolarization experiment while i45° polarization 
is the best choice when dual charnel measurements are to be made. 
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5.4 Polarization Diversity and Frequency Reuse 

It is the intent of this section to answer the questions relating 
to polarization effects which were posed in the Introduction. First, 
a polarization does exist for which the average attenuation during 
rainfall is a minimum. Figure 5.4-1 shows the attenuation on a 1.43 
km path at 19.3 GHz as a function of tilt angle. If it is assumed that 
positive and negative canting angles are equally likely then the time 
average canting angle is zero and the average tilt angle is equal to 
the polarization angle of the transmitted field. Then 0° and 90° tilt 
angles represent vertical and horizontal polarization respectively. 

From Figure 5.4-1, the attenuation is minimum for a tilt angle of 0°. 
Therefore, on the average, vertical polarization will suffer the least 
attenuation. 

What two orthogonal polarizations exhibit the least cross polar- 
ization interference? This question can be answered with reference to 
Figure 5.3-1 and with the assumption that the magnitude of the canting 
angle is about 15°. For a short path such as illustrated by the lower 
curve, vertical and horizontal polarizations produce the least cross 
polarization interference. For a canting angle of 15° the tilt angle 
for vertical polarization is 15° and the cross polarization level is 
-22 dB. The tilt angle for horizontal polarization is 90° - 15° or 75° 
and the corresponding level is -17 dB. For a 45° tilt angle the cross 
polarization level' is about -13 "dB. Clearly, horizontal and vertical 
polarization produce the least mutual interference for a 1.43 km path. 
However, this advantage decreases as the path length increases. For 
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longer paths, such as illustrated by the upper curve in Figure 5.3-1, 
the peak of the cross polarization versus tilt angle curve shifts toward 
horizontal polarization (90° tilt angle). This means that the cross 
polarization level for horizontal polarization is increasing more for 
a given increase in path length than is the level for other polarizations. 
This effect quickly eliminates the advantage of using vertical and 
horizontal polarization (at the same frequency) as two separate com- 
munications channels. 

This shift of the peak toward horizontal polarization is supported 
by Shimba [31] who made measurements of cross polarization interference 
using horizontal polarization at a frequency of 19.1 GHz, He observed 
a cross polarization level of -10 dB with a rain rate of 140 mm/hr on 
a 4.3 km path. Referring to the top curve in Figure 5.3-1, this implies 
a canting angle of about 8° since a -10 dB cross polarization level 
occurs for a tilt angle of 82°. A canting angle of 8° is within the 
expected range. 

The question of how long a path can be if polarization multi- 
plexing is to be used is answered in the next two figures for the two 
polarization configurations discussed here. Figure 5.4-2 shows the 
cross polarization level as a function of path length for a tilt angle 
of 60°. A tilt angle of 60° is the worst case condition for a +45° 
polarization multiplexed system and represents the 45° E field polar- 
ization plus a 15° canting angle. If a maximum acceptable cross polar- 
ization level of -10 dB is assumed, then the maximum path length is 
7 km for a 50 mm/hr rain rate. However, for a rain rate of 100 mm/hr 



Cross Polarization Level in dB 
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Figure 5,4-2, 


Cross i Polarization. vs. . Path. Length for a Tilt Angle 
of 60° and a Frequency of 19,3 GHz, 
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the maximum path length is reduced to 3 km and for 150 nim/hr it is only 
2 km. 

Figure 5.4-3 shows the cross polarization level versus path length 
for a tilt angle of 75°. This tilt angle is the worst case condition 
when a multiplexed system using vertical and horizontal polarization is 
employed and represents the 90° horizontal polarization minus a 15° 
canting angle. With the same -10 dB maximum acceptable cross polar- 
ization level, the maximum path lengths corresponding to rain rates 
of 50, 100, and 150 mm/hr are 11, 4.5, and 2.8 km respectively. Since 
rain rates of 150 mm/hr are not unconmon in many parts of the United 
States, the maximum path length would be limited to 2 km if t45° polar- 
ization multiplexing is to be employed and to 2.8 km if vertical and 
horizontal polarizations are used. 

The question as to the feasibility of using two polarizations 
which do not fade simultaneously to increase the reliability of a 
communications system is important from a practical standpoint. To 
answer this question one must realize that a system designed to minimize 
attenuation would of necessity employ vertical or near-vertical polar- 
ization since the attenuation is minimum for a tilt angle of 0°. For 
a constant rain rate the attenuation would vary as the canting angle 
varied from 0° to about 15°. From Figure 5.4-1 the difference in 
attenuation between a tilt angle of 0° and one of 15° is about 0.3 dB 
for the 1.43 km path or about 0.2 dB/km. Assuming vertical polarization, 
the minimum attenuation would occur for a canting angle of 0° and would 
be expected to increase by about 0.2 dB/km when the canting angle increased 
to 15°. For a 10 km path this represents an increase of, at most, 2 dB 
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above the minimum attenuation level. Hence, all that could be gained 
by employing another channel is 2 dB less attenuation for a 10 km path. 
This would indicate that polarization diversity is not feasible as a 
means of increasing resistance to fading in a communications system. 



SUMMARY 


Weather — especially precipitation — is a major consideration in the 
design of millimeter wave communications systems. The two major effects 
of weather are attenuation and depolarization of a transmitted wave due 
to rainfall. Most theoretical and experimental work to date has been 
concerned with the attenuation produced by rain while the depolarization 
effects have been comparatively ignored. The only theoretical model of 
depolarization prior to this work was the differential attenuation model 
proposed by Thomas. A few investigators had measured cross polarization 
levels during rain, but limitations inherent in their methods made an 
accurate comparison with theory virtually impossible. Hence, not only 
a new model of depolarizatiin by rainfall was needed, but also a new 
experiment conducted under conditions which would allow an accurate 
check of the theory. 

The work reported here satisfies both of these needs. A new 
theoretical model was developed based on the summation of the scattered 
fields from an ensemble of rain drops. The model forms the basis of 
a computer program which was written to calculate the expected attenuation 
and cross polarization levels during rainfall. Simultaneously, an 
experiment was begun at Virginia Polytechnic Institute and State Uni- 
versity to measure the effects of rain on millimeter waves. Data from 
this experiment are presented here and excellent agreement is shown 
between -the data and calculations based on the new model. 

Computations based on the new theory of depolarization developed 
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in this thesis allow several conclusions to be drawn regarding the 
influence of polarization on millimeter wave propagation through rain. 
The conclusions are: (1) the best polarizations to use for a depolar- 

ization experiment are i45° from the vertical, (2) vertical and hor- 
izontal polarizations should not be used for a depolarization experiment, 
(3) vertical polarization suffers the least average attenuation during 
rainfall, (4) Oguchi's attenuations and phase rotations for 19.3 GHz 
are correct, (5) the effective percentage of oblate drops assumed in 
an analysis is critical to the predicted cross polarization level, 

(6) polarization diversity is not feasible as a means of increasing 
resistance to rain-induced fading, and (7) the use of polarization 
multiplexing utilizing orthogonal polarizations is limited to very 
short path lengths. 
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APPENDIX 


NUMERICAL METHODS 


Before the scattering coefficients (S^’s) of Section 3.3.2 can 
be evaluated, the single drop scattering functions (f^'s) and the 
imaginary constant K must be determined* 

The f 's are computed using the scattering functions, and 
f^, which Oguchi has recently published [26], Assuming zero canting 
angle for a raindrop, f is the ratio of the scattered E field from 
the drop which is vertically polarized to the incident vertically 
polarized E field, and f^ is the ratio of the scattered E field from 
the drop which is horizontally polarized to the incident horizontally 
polarized E field. 

With reference to Figure A-l, E^ is the incident E field on the 
drop which has polarization 1 making an angle 0 with the vertical as 
shown. E* is decomposed into components E* and E* along the vertical 


= E* sin 0 


and horizontal axes respectively. 

E* = E ?■ cos 0 
v 1 

— ^ M g _g 

E^ and E^ are scattered by the drop producing E^ and E^. 

E 8 o f I 1 - f if cos 0 
V V v v 1 


K - - vi sln 9 


E^ and E^ have components E~ and E° in the polarization 1 direction 
which add to give the total scattered field which is 1 polarized. 


ss 


SS 


73 




75 


E® - E® cos 0 - f E* cos 2 0 

V^ V V 1 

E® = sin 0 = f h E^ sin 2 0 

E® = E® + E® = eJ (f cos 2 0 + f. sin 2 0) 

1 v x h x 1 ' v h 

— Q 

Then f^ is defined as the ratio of E^ to E^ 

S 1 2 2 
£ 11 " 51 * £ v cos 6 + f h sln 0 
E 1 

_g M g _g _ g 

E and have components E v and E^ in the polarization 2 direction 

2 2 

which add to give the total scattered field which is 2 polarized. 


E® = E® sin 0 = f E* sin 0 cos 0 
V2 v v 1 


E® = E® cos 0 «* f. E* sin 0 cos 0 
no h n 1 


e! = E S + E® = E* (f - f, ) sin 0 cos 0 
r 2 "2 


2 "v 0 ' “h 0 “1 v v h 

ss _ si 


Then f ^ is defined as the ratio of E^ to E^ 

S 2 

f 12 « «= (f v - f^) sin 0 cos 0 

S 2 

The coefficient f 22 i 8 obtained from the expression for by re- 
placing 0 with a where 

a - 0 - 90° 

l 2 2 2 
Then f„„ = — r *» f cos ct + f. sin a. 

II =i v h 

E 2 . - 


The determination of the constant K requires an evaluation of the 


I e 


-jkr 


summation 
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which is to be performed over all drops (in the plane at midpath) which 
are contained in the first Fresnel zone, and where r is the distance 
from the drop being considered to the receiver. With reference to 
Figure A-2, point 0 is the intercept of the line of sight with the 
plane of drops. Since only the relative phases of the waves from dif- 
ferent drops are significant, one may assume that the phase of the 
wave passing through point 0 is zero. The summation then becomes 


E e“ jk(r " R)2 

where the factor of 2 arises because the incident wave on the drop 
also travels the extra distance (r - R) . 

Then I e-J“< r * R > 2 = £ ^ 

If one divides the first Fresnel zone (radius = F^) into M intervals 
as follows 

r - F i 

p = /m - 1/2 — 0 < m < M 

” /H 


then the phase lag ip at any point is 


2kp ‘ , t F.‘ 

S_ = _ 2k _1_ ( _ 1/2 ) 

L L M v ’ 


From Section 3.2, 


- 2 xa(L-i) 

F i T 

XL/2(L-L/2) 

L 


X/2(L-L/2) 


at l - L/2 



plane of drops 
at midpath 


Transmitter 



Figure A-2. Geometry for Approximating Phase Lag of Scatter 


From Arbitrary Drop in the First Fresnel Zone 
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Then ¥ " - | <L - L/2) i (M - 1/2) 


5 < M - 1/2 > 


and the phase shift from point r to r ,, is constant at - — radians. 

m mfl M 

-jk(r - R)2 M -j ~ (m - 1/2) 

Now E e : ^ E e 

all M m=l 

drops 

since there are ^ drops in each interval where D is the total number 
of drops in the first Fresnel zone. 


M -j J (m - 1/2) 
But E e 


m=l 


+ *lk -Js(--i) 

e e E e 


= 21 V ~ J M P 

=* e l e 


-j 


JL 

2M 1 


p=0 


— j IT 
- e J 


-j — 
J 2M 


1 - e 
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“J H, 


"j ^ 

1 - e M 


" 1 ‘ -(-j *7 m) = " j 7 m for M large 

Then E e“ Jk(r " R)2 ; £(-j - M) - D 
all M TF 7T 

drops 


Hence, K ■ -J — 


The following computer program was written in FORTRAN IV incorporating 
the method developed in this thesis. The program is followed by a 
description of the variables and a discussion of convergence. 



COMPLEX FV, FH, FM, FX,MINUSJ , FXM, Sll , S22 , S12 , EM, EX, E 
REAL N 

COMPLEX CMP LX 

6 READ (5, 1)PATHL, RATE, THETA, N 

1 FORMAT (4F10.0) 

IF(N.EQ.O.)GO TO 5 
N MUST BE EVEN 
ALPHA=THETA- 9 0 , 

THETA=THETA/ 57.295 
ALPHA=ALPHA/57 . 295 
P=0.40 

DMODE=l .+0. 9*ALOG10 (RATE) 

FVR= (0. 1503+ ( (DMODE-2.) **1.256) *0.2226) *0.001 
FVI= (-0. 06914- ((DM0DE-2.)**1. 294 )*0, 19946) *0.001 
FHR= (0 . 167 7+ ( (DMODE-2 .) **1 . 306) *0 . 2986 ) *0 . 001 
FHI= (-0.07 75 7- ( (DMODE-2. ) **1.4 21)*0. 2653)*0. 001 
FV=CMPLX(FVR,FVI) 

FH=CMPLX(FHR,FHI) 

F11=FV*C0S (THETA) *C0S (THETA+FH*SIN (THETA) *SIN (THETA) 
F12= (FV-FH) *SIN (THETA) *COS (THETA) 

F22=FV*C0S (ALPHA) *COS(ALPHA)+FH*SIN (ALPHA)* SIN (ALPHA) 
R=PATHL/2. 

Rl=R— PATHL/ (2 . *N) 

R2=R+PATHL/(2.*N) 

V0LUME=3. 14159*0.0155* ( (R2*R2/2.-R2*R2*R2/ (3. *PATHL) ) 
1- (Rl*Rl/2.-Rl*Rl*Rl/ (3. *PATHL) ) ) 

TERMV=4.6*SQRT (DMODE) 

DR0PS=531 . *RATE/ (TERMV*DMODE*DMODE*DMODE ) *V0LUME 
MINUSJ=CMPLX(0. ,-l.) . 

S11=0.5*MINUSJ*F11*DROPS*2.*(2./3.14159)/R 
S22=0. 5*MINUSJ*F22*DROPS*2 . * (2 ./ 3 . 14159) /R 
S12=0 . 5*MINUSJ*F12*DR0PS*P*2 . * (2 . /3 . 14159) /R 
J=N 

EM=(1.+S11)**J 
E=CMPLX(0. ,0.) 

DO 2 1=1, J 

K=J-I 

L=I-1 

E=E+ ( 1 ,+Sll ) **K* (1 .+S 2 2 ) **L 

2 CONTINUE 
EX=E*S12 
EMR=REAL (EM) 

EMI=AIMAG (EM) 

EXR=>REAL(EX) 

EXI=AIMAG(EX) 

EMMAG=SQRT (EMR*EMR+EMI*EMI) 

EMPHAS=ATAN2 (EMI , EMR) 

EMPHAS=EMPHAS*57 . 295 
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EMDB=20.*AL0G10(EMMAG) 

EXMAG=SQRT (EXR*EXR+EXI*EXI) 

exphas=o. 

IF (EXMAG. EQ. 0. ) GO TO 9 
EXP HAS = AT AN 2 (EXI , EXR) 

EXPHAS=EXPHAS*57 .295 
EXDB=20 . *AL0G10 (EXMAG) 

EXPOL=EXDB-EMDB 
9 THETA=THETA*5 7.295 

WRITE ( 6 , 1 2 ) PATHL , RATE , THETA , N , P 

12 FORMAT(lHO, , PATHL= , ,F8.2,5X, , RATE= , ,F8.2,5X,'THETA= , ,F8.2,5X > , N= , > 
lFS^.SX/P^.FS^) 

WRITE (6,3) EMMAG , EMPHAS , EXMAG , EXP HAS 

3 FORMAT (1HO, 'EMMAG=' ,F10.5,5X, ' EMPHAS= * , F10 . 5 , 5X, ' EXMAG= ' ,F10.5,5X, 
l'EXPHAS=',F10.5) 

IF(EXMAG.EQ.O. )GO TO 10 
WRITE (6 , 4) EMDB , EXDB , EXPOL 

4 FORMAT (1HO, 'EMDB=' ,F10.5,5X, ’EXDB=' ,F10.5,5X, »EXPOL=' ,F10.5) 

GO TO 13 

10 WRITE(6,11)EMDB 

11 FORMAT (1HO , ’ EMDB= ’ , F10 . 5 ) 

13 WRITE (6, 8) 

WRITE (6, 8) 

WRITE (6, 8) 

8 FORMAT (1H0) 

GO TO 6 

5 STOP 
END 


Important variables in the above program are; 


PATHL - Path length in meters 
RATE - Rain rate in mm/hr 

THETA - Tilt angle of incident polarization in degrees 
N - Number of path segments (must be even) 

ALPHA - Tilt angle of cross polarization in degrees 
P - Fraction of effective oblate drops 

DMODE - Most frequently .occurring drop, size. (Bnpirical fit to Laws 
and Parson's data) 


FVR - Real part of single drop scattering function for vertical 
polarization (Empirical fit to Oguchi's data) 
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FVI - Imaginary part of above 

FHR - Real part of single drop scattering function for horizontal 
polarization (Empirical fit to Oguchi's data) 

FHI - Imaginary part of above 

FV - Oguchi’s scattering function for vertical polarization 

FH - Oguchi’s scattering function for horizontal polarization 

Fll - Scattering function, incident to incident polarization for 
one drop 

F12 - Scattering function, incident to cross polarization for one 

drop 

F22 - Scattering function, cross to cross polarization for one drop 

R - One-half the path length 

VOLUME - Volume of a segment of the path bounded by the first Fresnel 
zone 

TERMV - Terminal velocity of raindrops [8] 

DROPS - Number of drops of size DMODE in VOLUME 

511 - Scattering function, incident to incident polarization for 

one plane of drops 

512 - Scattering function, incident to cross polarization for one 

plane of drops 

S22 - Scattering function, cross to cross polarization for one plane 

of drops 

EM -Received signal, incident (main) polarization 

EX - Received signal, cross polarization 

EMMAG - Magnitude of received signal, incident polarization 
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EMPHAS 
EXMAG 
EXP HAS 
EMDB 
EXDB 
EXPOL 

The convergence of the above program as N is increased is shown 
in Figure A-3 where the cross polarization level is shown for a path 
length of 1.43 km, a rainrate of 150 mm/hr, a frequency of 19.3 GHz, 
and a tilt angle of 45°. As can be seen, convergence to within 1 dB 
of the final value is obtained for 1^20. 


- Phase of received signal, incident polarization 

- Magnitude of received signal, cross polarization 

- Phase of received signal, cross polarization 

- Received signal in dB, incident polarization 

- Received signal in dB, cross polarization 

- Cross polarization level in dB 




